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Abstract 

The delay time distribution (DTD) of type la supernovae (SNe la) from star formation is an important 
clue to reveal the still unknown progenitor system of SNe la. Here we report on a measurement of the 
SN la DTD in a delay time range of tja = 0.1-8.0 Gyr by using the faint variable objects detected in 
the Subaru/XMM-Newton Deep Survey (SXDS) down to i' ^ 25.5. We select 65 SN candidates showing 
significant spatial offset from nuclei of the host galaxies having old stellar population a,t z ^ 0.4-1.2, out 
of more than 1,000 SXDS variable objects. Although spectroscopic type classification is not available for 
these, we quantitatively demonstrate that more than ~80 % of these should be SNe la. The DTD is 
derived using the stellar age estimates of the old galaxies based on 9 band photometries from optical to 
mid-infrared wavelength. Combined with the observed SN la rate in elliptical galaxies at the local universe, 
the DTD in tia ^ 0.1-10 Gyr is well described by a featureless power-law as foitia) oc tf^ with a ~ —1. 
The derived DTD is in excellent agreement with the generic prediction of the double-degenerate scenario, 
giving a strong support to this scenario. In the single-degenerate (SD) scenario, although predictions by 
simple analytic formulations have broad DTD shapes that are similar to the observation, DTD shapes 
calculated by more detailed binary population synthesis tend to have strong peaks at characteristic time 
scales, which do not fit the observation. This result thus indicates either that the SD channel is not the 
major contributor to SNe la in old stellar population, or that improvement of binary population synthesis 
theory is required. Various sources of systematic uncertainties are examined and tested, but our main 
conclusions are not affected significantly 
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1. Introduction 

It is widely believed that type la supernovae (SNe 
la) are thermonuclear explosions of carbon-oxygen white 
dwarfs in binary systems, triggered when a white dwarf 
grows up to the Chandrasekhar mass by accretion from 
its companion (see Nomoto et al. 1997; Hillebrandt & 
Niemeyer 2000; Livio 2001 for reviews). However, the pro- 
genitor binary system leading to SNe la is still unknown, 
and there are two competing scenarios for the accretion 
process triggering SNe la. In the single-degenerate (SD) 
scenario, the accretion is from a non-degenerate compan- 
ion star (Whelan & Iben 1973; Nomoto 1982), while in the 
double-degenerate (DD) scenario, a merger of two white 
dwarfs results in a SN la (Iben & Tutukov 1984; Webbink 
1984). To reveal the progenitor is important not only for 
better understanding of this one of the brightest explo- 
sions in the universe, but also for controlling systematic 
uncertainties when SNe la are used as a standard candle 
to measure the expansion rate of the universe (Riess et 
al. 1998; Pcrlmutter et al. 1999). SNe la arc expected 
to have a wide range of delay time from star formation 
to supernova explosions, and the delay time distribution 



(DTD) can be used to discriminate the proposed progen- 
itor models, since different progenitor scenarios predict 
different DTDs. 

The DTD foitia) is equivalent to the SN la occurrence 
rate as a function of the la delay time tja, for a single- 
burst stellar population of a unit stellar mass. Hence, 
observational studies on SN la rates should be able to con- 
strain DTD (Mannucci et al. 2006 and references therein). 
One such approach is to examine the evolution of the cos- 
mic SN la rate density (Pain et al. 2002; Gal- Yam & 
Maoz 2004; Maoz & Gal- Yam 2004; Strolgcr et al. 2004; 
Dahlen et al. 2004; Barris & Tonry 2006; Forster et al. 
2006; Neill et al. 2006; Botticella et al. 2008; Oda et al. 
2008; Blanc & Greggio 2008). However, there is a degen- 
eracy between SN la DTD and the cosmic star formation 
history, and it is difficult to derive a strong constraint 
on DTD from the currently available data (Forster et al. 
2006; Botticella et al. 2008; Oda et al. 2008; Blanc & 
Greggio 2008). Another approach is to study the depen- 
dence of SN rate on the host galaxy properties, such as 
colors or spectral energy distribution (SED) (Mannucci 
et al. 2005; Scannapicco & Bildstcn 2005; Sullivan et 
al. 2006; Aubourg et al. 2007). This approach has al- 
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ready given some useful constraints on DTD, indicating 
the existence of SN la populations having both short 
0.1 Gyr) and long 10 Gyr) delay times. However, the 
functional shape oi foitia) has not yet quantitatively been 
constrained well. What has been done in previous studies 
is to assume simple mathematical functions for DTD or 
adopt theoretical DTD models to predict the distributions 
of some observational quantities such as host galaxy col- 
ors, and then compare them to the observed data. A clear 
next step is a measurement of DTD, i.e., to constrain the 
fimctional form of DTD directly from the observed data, 
rather than testing particular DTD functions or models. 

For a DTD measurement, one needs to reliably estimate 
the delay time of an observed SN la and the stellar mass 
of its host galaxy. The mean stellar age of the host galaxy 
is an indicator of the delay time, but it is unreliable when 
the stellar age distribution has a large dispersion, as in 
galaxies having extended star formation history. The ideal 
population for this purpose is old galaxies with an approx- 
imately uniform stellar age, which experienced major star 
formation episode in the past and have little or no star 
formation activity at the time of a SN explosion. Present- 
day elliptical galaxies are good examples, but they have a 
typical age of ^ 10 Gyr (Barber et al. 2007; Jimenez et al. 
2007), and hence we need to observe SNe at higher red- 
shifts to get a sample of SNe having shorter delay times. 
However, both detection of SNe and stellar age estimate 
of host galaxies become difficult at high redshifts. 

The Subaru/XMM-Newton Deep Survey (SXDS, 
Furusawa et al. 2008), which is the deepest survey among 
those wider than ~1 deg^ with a broad coverage of vari- 
ous wavelengths, provides a unique data set for this pur- 
pose. In optical bands, this field has been observed re- 
peatedly with some time intervals, and a systematic vari- 
able object search has been performed (Morokuma et al. 
2008a, b), leading to detection of more than 1,000 vari- 
able objects down to a limiting variability magnitude of 
nii' ~ 25.5 (AB). The majority of them are high-redshift 
supernovae and active galactic nuclei (AGNs). This data 
set includes many passively evolving old galaxies at red- 
shift z ~ 1, which are believed to be the direct ancestors of 
the present-day elliptical galaxies (Yamada et al. 2005). 
We can estimate the stellar age and stellar mass of host 
galaxies by photometric rcdshift calculations using the 
rich photometric data in the 9 band filters {BVRc i'z'JK, 
3.6, and A.bfim). The required accuracy for the age esti- 
mate in this work is about a factor of two, and we will 
show that such an accuracy can be achieved for the galax- 
ies used in this work, by a variety of tests for the age 
estimates. 

In this paper we measure the SN la DTD by select- 
ing the SXDS variable objects found in old or passively 
evolving galaxies at z ~ 0.4-1.2. Spectroscopic SN type 
confirmation is not available for the faint SXDS SN candi- 
dates. It is difficult to construct a complete high-redshift 
SN sample with spectroscopic type classification under ho- 
mogeneous conditions of spectroscopic observation, and 
contamination of SNe with no or poor spectroscopic data 
is one of the most challenging sources of uncertainty in 



SN rate studies (Strolger et al. 2004; Neill et al. 2006; 
Sullivan et al. 2006; Poznanski et al. 2007). Here, instead 
of using spectroscopic information, we select variable ob- 
jects in old galaxies with spatial offset from the galactic 
centers, and hence the majority of them are expected to be 
SNe la. In fact, we will quantitatively demonstrate that 
more than 80 % of them should be SNe la, based on the 
properties of the SN candidates and their host galaxies. 

It should be noted that selecting SNe only in a par- 
ticular type of galaxies docs not induce bias in the DTD 
estimates, since wc measure SN la rate normalized by stel- 
lar mass of host galaxies having the same type. This is in 
contrast to measurements of total cosmic SN rate density, 
in which selection of any particular galaxy type obviously 
leads to an underestimate of the total rate. On the other 
hand, we cannot measure DTD at short delay times ob- 
viously because of selecting old galaxies. We will present 
a measurement of DTD in a range of tia — 0.1-8.0 Gyr, 
corresponding to the distribution of mean stellar ages of 
the old galaxies used in this work. Our DTD measure- 
ment will be supplemented by the SN la rate measured 
for elliptical galaxies in the local universe, to obtain DTD 
in iia ^ 0.1-10 Gyr. The derived DTD will then be com- 
pared with a wide range of the existing theoretical DTD 
predictions, to get implications for the SN la progenitor. 

In section 2, we describe the SXDS data set and se- 
lection procedures of old galaxies and SN candidates. In 
section 3, we describe the formulations of DTD measure- 
ment and present the results. In section 4, we examine 
various systematic uncertainties in our DTD estimates. 
We then discuss about the implications for the SN la 
progenitor, from the comparison between the measured 
DTD and theoretical predictions (section 5). Summary 
and conclusions arc given in section 6 with some discus- 
sions. Throughout this paper we use the standard ACDM 
cosmological parameters of {h,il,M ,^a) = (0.7,0.27,0.73), 
where h = i7o/(100 km/s/Mpc), and magnitudes are given 
in the AB magnitude system unless otherwise stated. 

2. The SXDS SN la Candidates 

2.1. the SXDS Data and Photometric Redshift 
Calculations 

We utilize the rich photometric data set of SXDS in 
BVRc ^' z' bands obtained by the Subaru/Suprime-Cam 
(Furusawa et al. 2008), in JK bands obtained by the 
UKIDSS survey (Warren et al. 2007), and in 3.6 and 4.5 
fj,m bands obtained by the SWIRE survey (Lonsdale et 
al. 2004), with the Umiting magnitudes (3cr) of 28.4, 27.8, 
27.7, 27.7, 26.6, 24.1, 24.0, 23.1, and 22.4, respectively, 
for aperture diameters of 2 (optical and JK) and 3.8 (the 
Spitzer bands) arcsec. The point-spread-function size of 
the optical images is about 0.8 arcsec FWHM (1 pixel = 
0.2 arcsec). The total survey area used in the variable 
object survey (Morokuma et al. 2008a) is 0.918 deg^. 
Though the JK data is available only for 64.5% of the 
SXDS field, almost all (97.5%) of the SXDS field is covered 
by the Spitzer data. This is important for reliable stellar 
mass estimate, which is crucial for a study of SN la DTD. 
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Fig. 1. The properties of the old galaxies selected for the SN search. Upper left panel: the i'-band magnitude distributions of all 
the 16,492 old galaxies used in the analysis (solid), the 314 old galaxies with spectroscopic data (dotted), and the host galaxies of 
the 65 SN candidates (dashed). The dot-dashed line is the distribution of the i'-band variability magnitude of the 65 SN candidates. 
Upper right panel: the Ay distributions of the old galaxies (solid), non-old galaxies (dashed), and all (old plus non-old) galaxies 
(dot-dashed) with rrij/ < 24.0 and ms.s/jm < 22.8. Middle and lower panels: the distributions of (f*), tga/rgp, redshift, and stellar 
mass of all the old galaxies (solid) and the host galaxies of the SN candidates (dashed). 
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Therefore we select relatively bright 45,374 galaxies with 
niir < 24.0 and ma.g^m < 22.8, for reliable photo-z and 
stellar mass/age estimations. Stars have been removed by 
the two-color plot (to_r^ — mi' versus to_r^ — ms g^m) as in 
Morokuma et al. (2008a). 

Photometric redshift (zph) calculations are performed 
by using the publicly available code hyperz (Bolzonella, 
Miralles, & PeU 2000), with the GALAXEV hbrary 
(Bruzual & Chariot 2003) of stellar population synthe- 
sis models based on the Padova 1994 evolutionary tracks. 
In our baseline analysis, we use the library assuming the 
Salpcter initial mass function (IMF) in the mass range 
of 0.1-100 Mq and a supersolar abundance of Z = 0.05, 
since massive galaxies with old stellar population are ex- 
pected to have high metallicity as known for the local 
elliptical galaxies (Barber et al. 2007; Jimenez et al. 
2007). We use 7 models of star formation (SF) history 
having exponentially decaying star formation rate (SFR) 
with the exponential time scales of tsf = 0.1, 0.3, 0.7, 
1, 3, 5, and 15 Gyrs. We also add a model with a con- 
stant SFR, and hence 8 SF history templates are used 
in total. Attenuation of galaxy spectra by dust is taken 
into account with the Calzetti law (Calzetti et al. 2000) 
within a range of Ay — 0-2, where Ay is the restframe vi- 
sual band attenuation. The minimum photometric error 
is set to be ±0.05 mag, to include systematic errors for 
bright objects having negligibly small statistical photo- 
metric errors. Then the best-fit redshift, SF history tem- 
plate, galaxy age, and Ay are calculated for each galaxy. 
The age survey range is limited not to exceed the age of 
the universe at a given redshift. 

2.2. Old Galaxy Selection 

We select old or passively evolving galaxies by requir- 
ing tga/TgF > 2.3, where iga is the galaxy age estimated 
by the hyperz code, i.e., the time elapsed from the be- 
ginning of the template SF history. This value is chosen 
so that galaxies selected by this criterion have already 
formed more than [1 — exp(— 2.3)] = 90% of stars that the 
galaxies would form in all the history. We also set a con- 
straint of Ay < 1.0 to avoid dusty objects. After selecting 
galaxies within 0.4 < Zph < 1.2, which is a typical redshift 
range expected for SNe la detectable in SXDS, there re- 
main 16,492 galaxies. We call these galaxies simply "the 
old galaxies" in this work for convenience. 

Stellar masses are calculated for each of the old galax- 
ies from the results of the photo-z fit. In this work, we 
define the stellar mass of a galaxy, M*, as the integra- 
tion of SFR up to its age. Because dying stars return a 
part of their mass into interstellar medium, this is not ex- 
actly the same as the mass locked up in stars at a given 
age, which changes with time even in a passive evolution 
phase. [The difference between the two is typically ^ 30%; 
see, e.g.. Fig. 7 of Sullivan et al. (2006).] It should be 
noted that the SFR-integrated mass is more appropriate 
for our purpose, because DTD should be normalized by 
the amount of star formation for a single starburst stellar 
population, which is constant against age. However, stel- 
lar mass estimations for observed galaxies are generally 



model-dependent. Therefore, when we derive DTD, we 
convert the stellar mass of our definition into Lk.o, which 
is the restframe K band luminosity at the age of 11 Gyr 
(a typical age of local elliptical galaxies; see e.g., Baber 
et al. 2007; Jimenez et al. 2007) after passive evolution, 
since restframe K band luminosity at a fixed large age is 
a good indicator of stellar mass. Then DTD will be given 
per unit Lk,o, and hence the difference in mass definition 
does not affect our DTD measurements. By this normal- 
ization, we can avoid some observational uncertainties in 
stellar mass estimates such as IMF (see also section 4.3), 
and a comparison between our results and other SN rate 
studies also becomes easier. 

The distributions of i' magnitude, Ay, mass- weighted 
mean stellar age (i,), iga/TsF, redshift, and stellar mass 
of these old galaxies are shown in Fig. 1. As expected for 
old galaxies, the distribution of Ay is peaked at Ay = 0, 
with a mean value of {Ay) = 0.31. The concentration to 
small Ay values is obvious when it is compared with the 
distribution of the non-old galaxies. To test the reliabil- 
ity of photo-z, we compare the results with the spectro- 
scopic redshifts (zgp) of the 314 old galaxies having ob- 
served spectra in Fig. 2. The agreement is within ±20% 
for the majority of the old galaxies. As a measure of de- 
viation that is not sensitive to outliers, we computed the 
median of 1.48| Az|/(1 + z), which is the same as the stan- 
dard deviation of Az/(1 + z) in the case of the Gaussian 
distribution. The result is 0.035, and this is comparable 
with other photometric redshift studies at similar redshifts 
(e.g., Ilbert et al. 2006). The magnitude distribution of 
the spectroscopic sample peaks at m^/ 22-23 and extends 
down to mil ^ 24 (Fig. 1), and hence our photometric old 
galaxy sample is restricted to the magnitude range where 
the spectroscopic calibration is possible. In the DTD anal- 
ysis below, we use spectroscopic redshifts when available, 
and otherwise use photometric redshifts. 

For visual demonstrations, we randomly selected 8 ob- 
jects as examples from our final sample of SN candidates. 
The properties of the 8 objects are summarized in Table 1. 
The host galaxy images and the SED fits of the 8 objects 
are shown in Figs. 3 and 4, respectively. The photometric 
errors are very small especially in the optical bands since 
we selected bright objects, but still the agreement between 
the template SEDs and the observed data is quite good, 
which is an encouraging result about the reliability of age 
and stellar mass estimates. 

2.3. Selection of the SN Candidates 

We now search for variable objects associated with these 
old galaxies. There are 1,040 variable objects selected 
by i' band variability in the catalog of Morokuma et al. 
(2008a). The SXDS field was monitored 8-10 times during 
about four years, with a typical time interval from a few 
days to one month within a year. The variable sources are 
selected if it is detected in the subtracted images of any 
possible pair of two different epochs. The detection effi- 
ciency of variable sources has been carefully determined 
by simulations. "The variability fiux" of a SN candidate in 
this work refers to the maximum of differential fluxes mea- 
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Fig. 2. Spectral redshift versus photometric redshift of the 
old galaxies selected for the SN la search. The dashed lines 
indicate the ±20% accuracy region. 

sured in time intervals corresponding to all possible pairs 
of the SXDS observing epochs. "The variability magni- 
tude" is the magnitude corresponding to this variability 
flux. 

We select variable objects within the detection isophotc 
(i.e., the area where the surface brightness is higher than 
the threshold level of the source detection) of the old 
galaxies, and showing significant offsets (> f.9 pixel = 
0.38") from the nuclei of galaxies. The mean area of de- 
tection isophote of the old galaxies is 409 pixels. Here, 
the nuclei of galaxies are simply defined by their sur- 
face brightness peaks. According to the simulations us- 
ing artificial objects (Morokuma et al. 2008a), this is a 
conservative limit to remove variable objects at nuclei of 
the host galaxies, i.e., AGNs. Morokuma et al. classi- 
fied the nuclear variable sources into SNe or AGNs based 
on the light-curve shape, and they found that such clas- 
sifications are broadly consistent with the X-ray informa- 
tion. However, the classification is rather uncertain espe- 
cially for the faintest variable objects near the detection 
limit. Therefore we conservatively reject all nuclear vari- 
able sources. We do not use the X-ray information in the 
process of SN candidate selection, since a small part (~ 
10%) of the survey area is not covered by the X-ray ob- 
servation. Galaxies detected in the X-ray band are only 
68 out of the 14,909 old galaxies observed in X-ray, which 
is a negligible fraction in our analysis. In fact, wc con- 
firmed that no object is detected in X-ray in our final SN 
candidate sample. 

We thus found 67 variable objects associated with the 
old galaxies. We examined the optical variability infor- 
mation of these objects by using the full Suprime-Cam 
data set of the SXDS. Wc found that, as expected, most 
of them are consistent with SN-like variability. Only 2 
objects show clearly AGN-like variability, i.e., flaring up 
more than twice during the four year observation dura- 
tion. These two are removed, leaving the 65 SN candi- 
dates. Such a contamination of AGNs is possible by a 
chance superposition of a normal galaxy and an unrelated 



AGN along the line of sight (e.g., the case of an AGN 
initially classified as SN 1999de, Gal- Yam et al. 2008). 
We found that two such events are reasonable from the 
number of AGNs in SXDS and the area covered by the 
old galaxies selected here (2.3% of the total survey area). 
We cannot exclude a possibility that a few more chance- 
superposition AGNs having SN-like light curves might be 
included in the 65 SN candidates. To check this, we ex- 
amine the offset distribution of these objects with respect 
to the surface brightness profile of galaxies as follows. 

2.4- Radial Distribution of the SN Candidates in Host 
Galaxies 

We performed ellipsoidal fits for each galaxy with the 
Sersic type radial profile (Sersic 1968). The effective ra- 
dius along the major and minor axes, and the Sersic index 
are thus derived, taking into account the seeing. We then 
calculated the fraction, //, of z'-band light enclosed within 
the ellipsoidal radius to a variable object, compared with 
the total flux within the detection isophote (i.e., isophotal 
flux). Here, the nuclear region within the offset threshold 
(1.9 pixel radius) is excluded in the // calculation. If we 
are selecting SNe la physically associated with the galax- 
ies, we expect that the distribution of // is roughly uniform 
between 0-1, since SNe la are expected to trace the host 
galactic light approximately (Kelly et al. 2007; Forster 
& Schawinski 2008). The distributions of the offsets from 
nuclei and /; are shown for the 65 SN candidates in Fig. 
5, and in fact we find an almost uniform distribution for 
fi- 

There is a deficit at // ^ 0.8, and it indicates that our 
profile fitting is not perfect or supernova locations do not 
exactly trace galactic j'-band light. However, it is clear 
that the SN candidates trace the galactic light reason- 
ably well; and hence the majority of these must be SNe 
physically associated with host galaxies, rather than the 
chance superpositions of background AGNs. It should be 
noted that, if there is a significant contamination of AGN 
chance superpositions, we expect a distribution biased to 
larger values of /;, which is opposite from the observed 
trend. We thus define these 65 objects as the final SN 
candidate sample, and the distributions of variability and 
host galaxy magnitudes are shown in Fig. 1. The images 
of the randomly selected eight examples are shown in Fig. 
3. 

In the final SN sample, 7 SN candidates have measured 
spectroscopic redshifts of host galaxies. We confirmed 
that the difference between Zgp and ^ph is within ~20 % 
for all of these, as expected from Fig. 2. The fraction 
of available spectroscopic redshifts, 7/65, is considerably 
higher than 314/16492 for the old galaxies. This is mainly 
because the host galaxies of the SN candidates are system- 
atically brighter than the typical old galaxies (see Fig. 1), 
reflecting the fact that SN rate is roughly proportional 
to host galaxy luminosities. Another possible reason is 
that some variable objects were followed up by spectro- 
scopic observations as a part of the Supernova Cosmology 
Project in collaboration with the SXDS project. 
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Table 1. Properties of the eight examples randomly selected from the final 65 SN sample. The redshifts are spectroscopic when 
available (objects 2 and 4), and otherwise photometric. Magnitudes of host galaxies and SN variability are in i' band, and galaxy 
age (iga), the exponential decay time scale of SFR (rgp), and mass-weighted mean stellar age ((t*)) are all in Gyrs. 
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Fig. 3. The images of the eight examples randomly selected from the final SN candidates. The left panels are the 3-color composite 
images of the host galaxies using BVi' bands. The middle panels are color contours of the host galaxy i' band flux at the time of 
the maximum luminosity of the SN candidates, and the right panels are those of subtracted images showing the variability of the 
SN candidates. All panels are centered at the surface brightness peak of the host galaxies, and the locations of the SN candidates 
are indicated by white crosses. The properties of the eight objects are summarized in Table 1. 
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Fig. 4. The results of the SED fits for the old galaxies by the photometric redshift calculations, for the 8 example objects listed 
in Table 1. The SED is given by energy flux per unit wavelength in arbitrary units. The open squares are the observed flux. The 
horizontal error bars indicate the approximate width of band filters, and the vertical error bars are la flux errors (difficult to see in 
optical bands). The flux errors include only statistical errors, while the minimum error is set to be ±0.05 mag in the photometric 
redshift calculations to take into account the systematic errors. 
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Fig. 5. Loft panel: the distribution of the spatial offsets of the 65 SN candidates from the nuclei of their host galaxies. Right panel: 
the distribution of the fraction of j'-band galactic light /; enclosed within the ellipsoidal radius to a SN candidate compared with 
the isophotal flux of its host galaxy. 



3. Measurement of Delay Time Distribution 

3.1. The Basic Formulations 

Now we estimate the DTD from the 65 SN candidates. 
Later (section 4) we wih demonstrate that the majority of 
them must be SNe la rather than core-collapse (CC) SNe, 
but here we first estimate the DTD, taking into account 
the contamination of CC SNe. We implicitly assume that 
DTD is universal for all stellar populations, as often as- 
sumed in studies on the SN la DTD. It is possible that 
DTD is dependent on the properties of stellar population 
such as metallicity (e.g., Kobayashi et al. 1998), and in 
such a case our measurement of DTD applies only to the 
old galaxies selected here. This point should be kept in 
mind when a comparison is made with theoretical predic- 
tions. 

Since we know the star formation history of the tem- 
plate galaxy evolution model in the hyperz fit, we can 
estimate the mass-weighted mean stellar age of galax- 
ies, (t*). Since we selected old galaxies, in which most 
of stars have formed in a starburst in the past, (i*) is 
expected to be a good estimator of the la delay time, 
^la. (By the same reason, the difference between the 
mass- weighted and luminosity-weighted mean stellar ages 
should be small.) The DTD function foitia) (per unit 
stellar mass and unit delay time) can be estimated as 
the SN la rate per unit stellar mass in these galaxies. 
Therefore the first rough estimate of the mean DTD func- 
tion in a given bin of ti < < can be obtained by 
solving the following equation: 

.^Ia,exp + .^CC.cxp = .^obs , (1) 

where A'la^cxp and iVcc.cxp are the expected numbers of la 
and CC SNe, respectively, and A'obs is the number of the 
observed SN candidates associated with the old galaxies 
satisfying ti < (t^,) < tu. We can calculate A'ia,cxp as: 

./Vla.exp = .foy^/off,! Mt,i ^'^^^ , (2) 



where M^^^ and Zi are the stellar mass and redshift of i-ih 
galaxy, fos^i the fraction of light (in i' band) outside the 
threshold offset (1.9 pixel radius), and the summation is 
for all the old galaxies satisfying ti < (i*) < In this 
work, we set the shortest delay time bin to be 0.1-0.25 
Gyr. This is reasonable because the minimum value of 
(t*) possible for the old galaxies selected by our criteria is 
0.16 Gyr corresponding to the case of tga = 0.23 Gyr and 
'''SF =0.1 Gyr, with an age dispersion of about 0.1 Gyr. 
On the other hand, the number of the old galaxies with 
(t*) > 8.0 Gyr is very small because they are located at 
z > 0.4, and hence we set the upper boundary of our DTD 
measurement to be 8.0 Gyr. 

The visibility time Tv (also often called as the control 
time), which is the total integrated time during which a 
SN can be detected by the SXDS observations, can be 
calculated as: 

/oo 
Pdot{t,z)dt , (3) 
-oo 

where t is the explosion time of a supernova relative to the 
SXDS observation campaign (in observer's frame), and 
P(iQtit,z) is the detection probability of the supernova at 
redshift z. The factor (1 -I- z,)"^ in eq. (2) corrects the cos- 
mological time dilation effect. This is the standard quan- 
tity in SN rate analyses, and it can be calculated if the SN 
light curve is given in the observing band filter. We ap- 
ply the standard light curve and color evolution of SNe la 
used in Oda & Totani (2005) , taking into account the dis- 
persion of the peak B band luminosity and the light-curve 
versus peak luminosity relation. We assume an extinction 
corresponding to the reddening E{B — V) = 0.05, which is 
reasonable for old galaxies and often used in rate studies 
of SNe la [see Oda et al. (2008) and references therein]. 
We assume the standard Milky- Way extinction curve. It 
should be noted that the Calzetti attenuation law used in 
the photo-z calculations is for the synthesized flux from 
a galaxy, which is different from the extinction law for a 
single star in a galaxy. Therefore Ay of a galaxy esti- 




Fig. 6. The visibility time Ty (in observer's frame) of SNe 
in SXDS as a function of redshift. The solid curve is for the 
standard case of SNe la, while the dot-dashed curve is for CC 
SNe. The dashed curves are for SNe la, but the luminosity is 
shifted by ± 0.4 mag. The visibility time is slightly different 
in the different subfields of SXDS, and those in the SXDS-C 
field (Morokuma et al. 2008a) are shown here. 

mated by photo-z calculation cannot be directly related 
to extinction of a SN in the galaxy. The effect of chang- 
ing extinction of SN flux will be examined in section 4. 
Finally, the exact procedure and efficiency of variable ob- 
ject detection in SXDS (Morokuma et al. 2008a) are also 
taken into account to calculate Pdet- The calculated vis- 
ibility time is shown in Fig. 6. Note that this visibility 
time includes the multiplicity of many SXDS observations 
during four years, and not for a single epoch observation. 

The expected number of CC SNe, iVcc,oxp, can be cal- 
culated as 

Ncc.cxp fos.i -01 fee '^^''^^^^'■^ , (4) 

where i]} is SFR (mass per unit time) of a galaxy, and 
/cc the production efficiency of CC SNe per unit mass 
of star formation. The SFR ?/» is estimated by the ob- 
served B band flux corresponding to the rcstframc UV 
luminosity, with the conversion factor calculated by the 
stellar population synthesis model used in this work. The 
factor /cc is calculated by assuming that all stars heav- 
ier than 8 Mq produce CC SNe with the assumed IMF. 
The visibility time TV,cc is calculated with a standard 
mixture of light curves of various types of CC SNe, as in 
Oda & Totani (2005). We assume E[B -V) = 0.15 for 
CC SNe, which is a typical value for CC SNe found in 
nearby galaxies (Oda ct al. 2008 and references therein). 
This value is larger than that assumed for SNe la, but it 
is reasonable since CC SNe are expected to occur in star 
forming regions where dust abundance is generally high. 
Since we selected old galaxies, we found A^cc.exp = 11-4 
(in all the considered range of = 0.1-8 Gyr) and this 
is small compared with A'obs = 65. (See also Table 2 for 
numbers in each bin.) 

The DTD estimated by this simple formulation is shown 
by open squares in Fig. 7. The error bars are statistical Icr 



Fig. 7. The SNe la DTD, foitla.), pcr unit delay time ti^ 
[century"^] for a single starburst population whose total 
iC-band luminosity is 10^^ Lfc q at the age of 11 Gyr. The 
filled squares are the final observational estimates by this work 
based on the baseline analysis, and the error bars are statis- 
tical la errors. The open squares are the same but using a 
simpler method to estimate the delay time. The time bins are 
the same as those for the filled squares, but the open squares 
are shifted in time by +0.03 dcx in this plot for presenta- 
tion. The open circle is DTD inferred from the SN la rate in 
elliptical galaxies in the local universe (Mannucci et al. 2005). 

errors, which are calculated by the confidence limits of the 
small number Poisson statistics (Gchrcls 1986). As men- 
tioned in section 2.2, the normalization by stellar mass 
has been converted into that by the K-hand luminosity 
Lx,o, i.e., the luminosity at the age of 11 Gyr after pas- 
sive evolution. The conversion is performed by calculat- 
ing M^/Lk,o from the SED templates used in the photo-z 
calculations. For the template of rgp =0.1 Gyr, we find 
M*/Lk.o = 1-8 [Mq/Lk,q]. Here, we used {V - K)q = 
+ 1.49 (Cox 2000) and fx = 4.17 x 10"" erg ciR-^s-^k-^ 
for i^T = (Zombeck 2007) in the Vega magnitude system. 

We assume that the stellar age of nearby elliptical galax- 
ies is 11 Gyr, based on the estimates for the SDSS early 
type galaxies (Barber et al. 2007; Jimenez ct al. 2007, and 
see section 4.5 for more details). Hence, the observed SN 
la rate per unit K luminosity in nearby elliptical galaxies, 
0.035l^;5i?(/i/0.75)2 century-i {IO^^Lk^q)-^ (Mannucci 
et al. 2005), gives an estimate of /_d(11 Gyr), which can 
directly be compared with our estimates. This value is 
based on the rate measurement of CappcUaro et al. (1999) 
assuming Hq = 75 km/s/Mpc, and it is translated into a 
value for ft. = 0.7 used in this work. 

Though these DTD results are derived by simply as- 
suming iia = (t* ) , they are not significantly different from 
our final results shown by the filled squares (sec below), 
and well described by a power-law, foitia) oc tia" with 
a ~ —1 in 0.1-11 Gyr. However, to get a more accurate 
estimate, we make a correction as follows. 

3.2. Correction for the Delay Time Estimates 

In the above formulation, we simply estimated tia of an 
observed SN la by (t*) of its host galaxy. However, in 
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Table 2. The DTD measurements by the basehne analysis. The DTD results are shown in [century"^] for a single starburst 
population whose total K band luminosity is lO^'^Lj^ q at an age of 11 Gyr. The errors are statistical la. The number of the old 
galaxies (A^gai); that of the detected SN candidates (A^obs)i the expected numbers of the prompt SNo la (A''pja_cxp) a^nd CC SNe 
(-^CC.exp) are shown. The last column shows the mean of atj^ of galaxies in a time bin, where crtj^ is the standard deviation of the 
probability distribution of tia in a galaxy. 



fia bin [Gyr] 


DTD /z,(tla) 






^ pla,cxp 


^CC.cxp 


{^tj [Gyr] 


0.1-0.25 


2 89+^-''" 
z.oy_]^ 20 


3719 


12 


1.0 


2.7 


0.05 


0.25-0.5 




3086 


14 


0.92 


3.1 


0.14 


0.5-1.0 


0.8410-3° 


4260 


19 


0.88 


3.0 


0.30 


1.0-2.0 


0.49l°j^ 


4764 


19 


0.66 


2.3 


0.50 


2.0-4.0 




564 


1 


0.051 


0.20 


1.2 


4.0-8.0 


00+1°® 


99 





0.007 


0.032 


2.1 



reality, there is a probability distribution of <ia, which is 
determined by the DTD and star formation history il}{t), 
where t = tga — tja is the time elapsed from the beginning 
of star formation in the host galaxy. The expectation 
value of the la delay time, (tia), is exactly the same as (t*) 
only when foiiia) is constant against tja. The expectation 
value is generally given by: 

Therefore, a more accurate estimate is obtained by the 
same formulation but with a summation over galaxies sat- 
isfying ti < (tia) < t^. However, this integration diverges 
when tia if a < — 1. The possibility of a significant 
population of the prompt SNe la with tjg, ^ 0.1 Gyr has 
been discussed in recent years (Scannapieco & Bildsten 
2005; Mannucci et al. 2006; Sullivan et al. 2006; Aubourg 
et al. 2007). Such a prompt population may also affect 
the estimate of the delay time. Therefore, we separate 
the "prompt" population (defined by tja < = 0.1 Gyr) 
from the "delayed" ("tardy") population (tia>tp). Then, 
the DTD of the delayed component can be estimated by 
a modified version of cq. (1), 

-^Ia,Gxp + -^pla.oxp + ^CC.cxp = ^obs , (6) 

where iVia,oxp is now for the delayed population and the 
integration of eq. (5) is from tia = tp to tga,. The expected 
number of the prompt la, A^pia,exp, can be calculated from 
SFR in a similar way to get A^cc.cxp, if the integrated DTD 
at < tia < tp is given. 

An obvious difficulty in this new formulation is that 
we need to know the DTD itself beforehand to make a 
measurement of DTD, because (tia) and A^pia.cxp depend 
on /_D(tia). In fact, this circularity problem can be solved 
easily; first we estimate fo with an initial guess of foitin), 
and then we can iterationally repeat this process with the 
new estimate of foitia.)- The first DTD estimate with 
a constant fu prior (open squares in Fig. 7) is well de- 
scribed by a power-law (oc iia") at tia > tp, and hence we 
assume this form of fo for the delayed population at the 
later iterations. For the prompt population, we assume 
that fo is constant at ija < tp and it is continuously con- 
nected to the delayed component at <ia = tp, for every 
iteration. This is reasonable from the DTD predictions 



by the stellar evolution theory (section 5), and the effect 
of this assumption on the DTD measurement will be ex- 
amined in section 4. 

We get sufficiently convergent results by just a few it- 
erations of this procedure, and they are shown by filled 
squares in Fig. 7 as our final results. Since we selected 
old galaxies, (t*) is already a good estimate of (iia), and 
this is why the correction is not large. The expected num- 
ber of the prompt population, iVpia.cxp, is 3.6 for the en- 
tire range of tia = 0.1-8 Gyr, which is just 5.5% of the 
65 SN candidates. The derived DTD, as well as im- 
portant quantities such as A^pia,exp, A^cc.cxp, and iVobs, 
are summarized in Table 2 for the six time bins in tia 
= 0.1-8 Gyr. The best fit power-law to the measured 
DTD is foil Gyr) = 0.55t°;i? century-i(lOi°Lx,o,0)-' 
and a = — 1.08+o;^5. Here, we included the data point at 
tia = 11 Gyr in the fit. 

4. Examination of Systematic Uncertainties 

4-1. Are the SN candidates really SNe la? 

In our statistical estimation of DTD, it is not necessary 
to prove that all of the SN candidates are SNe la. Since 
we have already shown that the selected variable objects 
are tracing the light profile of host galaxies, the major- 
ity of them must be SNe la or CC SNe. What we need 
to demonstrate is, then, that about 82% of the 65 candi- 
dates are actually SNe la, as inferred from the estimate of 
A^cccxp = 11.4. 

We first examine the i' band variability magnitude 
'7ivar,i' of the 65 SN candidates, and compare with the 
brightest magnitude TOmin,!' that is possible for SNe la 
corresponding to the peak luminosity. We calculated this 
magnitude for each SN candidate assuming the standard 
type la light curve with the mean peak B magnitude used 
in Oda & Totani (2005). We find that most (61/65) of the 
SN candidates have mvar.i' > Tninins' as expected, and all 
objects satisfy TOvar,i' > ™min,i' — 0.53. The four objects 
with mysx_ii < mmin,i' are reasonable considering the dis- 
persion of the peak B magnitude of SNe la (~0.4 mag). 
On the other hand, about half of the SN candidates have 
the variability magnitudes brighter than TOmin,i' + 1, and 
such variability cannot be explained by CC SNe, because 
most of them are fainter than SNe la by 1-2 magnitudes 
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Table 3. The results of the Kohnogorov-Smirnov tests for 
the redshift and host galaxy stellar mass distributions of the 
SN candidates. The chance probabilities of getting the ob- 
served distributions are shown. (See Fig. 8 for the graphical 
presentations of these distributions.) The columns 2-4 show 
the results when all SNe are assumed to be delayed SNe la, 
prompt SNe la, and CC SNe, respectively. The last column 
shows the results when the three populations are mixed with 
the expected numbers (45.2, 3.6, and 11.4) in our baseline 
analysis. 





delayed la 


prompt la 


CC 


mix 


redshift 


0.07 


0.81 


0.05 


0.40 


stellar mass 


0.05 


0.07 


0.04 


0.48 



(see e.g., Oda & Totani 2005). 

We already made the DTD estimates, and we can 
predict the expected distribution of redshifts and stel- 
lar masses of host galaxies for the SN candidates, be- 
cause we can calculate the expected number of SNe la 
in each galaxy. Comparison of these expected distribu- 
tions with those observed provides us with an important 
consistency check of our DTD estimates. Figure 8 shows 
such comparisons, and the predictions based on our final 
DTD estimates are in nice agreement with the data. In 
the calculation of the expected redshift and stellar mass 
distributions, the light curve information of SNe is in- 
cluded through the visibility time. Therefore, we do not 
expect such an agreement, if our estimate of la/CC ratio 
is wrong or there is a significant contamination from any 
non-SN objects. In fact, the predicted distributions are 
in serious contradiction with the data if we assume that 
all the SN candidates are CC SNe. The agreement be- 
tween the expected and observed distributions are quan- 
titatively tested by the Kolmogorov-Smirnov test, and the 
results are given in Table 3. Acceptable fits are obtained 
both for the redshift and stellar mass distributions only 
when we make an appropriate mix of the delayed/prompt 
la and CC SNe with the relative proportions predicted by 
our DTD estimates: iVexp = 45.2 (delayed la), 3.6 (prompt 
la), and 11.4 (CC). These results then give a strong sup- 
port to the reliability of our DTD estimates. 

It should also be noted that we selected all variable ob- 
jects associated with the old galaxies with a significant 
offset from the nuclei. Only 2 objects were removed by 
AGN-like variability, which is a negligible number com- 
pared with the final 65 candidates. Therefore, our DTD 
measurement gives a conservative upper limits to the true 
DTD, and the derived DTD would always be overesti- 
mates if there was any unknown contamination from non- 
SN objects. 

Finally, we repeated our DTD estimates with a more 
stringent criteria for the old galaxies, iga/TsF > 3.0, to 
select more efficiently SNe la rather than CC SNe. In 
this case, the number of the SN candidates is reduced to 
44, but now A^cc.cxp = 5.4 and A^pia,oxp = 1-6, and hence 
the expected la (dclayed-|-prompt) fraction is increased 
to 88% from 82% in the baseline analysis. The DTD es- 
timates and the power-law fit in this case are shown in 



Table 4 and Fig. 9, and the difference from the baseline 
results is within the statistical uncertainties. It should be 
noted that the number of the SN candidates in the short- 
est bin (0.1-0.25 Gyr) is now only 5, because the possible 
minimum of mean stellar age is increased to 0.22 Gyr from 
0.16 Gyr in the baseline analysis. Therefore, although the 
DTD value in the shortest time bin shows a larger devi- 
ation from the baseline analysis than in other time bins, 
the statistical uncertainty is quite large. 

4-2. Estimates of (tja) and contribution from the prompt 
la population 

In our analysis, the delay time iia of each SN candi- 
date has been estimated simply by the expectation value 
of (tia) from the star formation history of the host galaxy. 
We have already shown that the difference between (tis) 
and the mean stellar age (t*) is small and does not seri- 
ously affect the final results. However, the probability dis- 
tribution of tia should have some dispersion around (tia), 
and if this dispersion is larger than the bin width of the 
DTD estimate, it could affect the results. The dispersion 
around (iia), ftia = ((^L) ~ (^la)^)^^^ is therefore calcu- 
lated for each host galaxy in a similar way to eq. (5), and 
their average in a tja bin is shown in Table 2. The disper- 
sions are typically about half of the bin widths, and hence 
the binning size of our DTD estimates is appropriate. 

In our baseline analysis, the amount of prompt SNe la is 
calculated by assuming a constant fnitin) at tia < tp, and 
in this case the fraction of the prompt la is about 20% 
of all SNe la when the DTD is integrated over a range 
of iia = 0-11 Gyr. However, since we selected old galax- 
ies, iVpia,oxp = 3.6 is just 5.5% of our 65 SN candidates. 
The expected fraction of prompt SNe la is the largest in 
the shortest delay time bin (0.1-0.25 Gyr), but it is still 
less than 10% (1.0/12, see Table 2). Therefore, our re- 
sult will not be seriously affected by changing the prompt 
fraction, unless we assume an extremely higher prompt la 
rate (e.g., by a factor of about 10) than that assumed in 
the baseline analysis. We show the results of the DTD 
estimate when the prompt la population is enhanced by 
a factor of 2.5 from our baseline analysis in Table 4 and 
Fig. 9. In this case, the prompt fraction in all SNe la 
integrated over tia =0-11 Gyr becomes ^ 50%, which is 
a typical value discussed in recent papers about prompt 
SNe la (e.g., Mannucci et al. 2006; Sullivan et al. 2006). 
Our DTD estimate for the delayed SNe la is not seriously 
affected. 

4-3. Systematic Uncertainties in Stellar Age and Mass 

The estimates of stellar mass and age of the host galax- 
ies are crucial in our DTD estimates. Since we selected old 
or passively evolving galaxies in a relatively bright mag- 
nitude range of m^/ < 24, the estimates are expected to be 
easier than general studies of high-z galaxies. However, 
these estimates may sensitively depend on e.g. the SF 
history, IMF, and metallicity assumed in the photometric 
redshift calculation. Therefore the systematic uncertain- 
ties must be carefully examined. 

We test several age and mass estimations with different 
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Fig. 8. The cumulative distributions of redshift and stellar mass of host galaxies for the SN candidates. The thick solid lines are 
the observed distributions of the 65 SN candidates. The dashed and dot-dashed curves are the expected distributions for the delayed 
(tia > 0.1 Gyr) and prompt (tja < 0.1 Gyr) components of SNe la, respectively, and the dotted curve is that for CC SNe. The 
thin solid curve is our best guess, which is the sum of the three components with the expected numbers in the baseline analysis: 
A^exp =45.2 (delayed la), 3.6 (prompt la), and 11.4 (CC). The results of the Kolmogorov-Smirnov test for these distributions are 
given in Table 3. 



Table 4. Examination of systematic errors in the DTD measurements. Various DTD results are shown when the analysis method 
is changed from the baseline analysis. The last two columns show the best-fit power-law DTD, foiila,) = /D,lGyr(tla/l Gyr)". 
The "tia = (t*)" results are obtained by estimating tja simply by mean stellar age (t*) of the host galaxies (open squares in Fig. 
7). The "tga/Tgp > 3.0" results are obtained when a more strict criterion of the old galaxies is applied than the baseline analysis. 
The "prompt x2.5" results are obtained when the fraction of the prompt la population is increased by a factor of 2.5 from the 
baseline analysis. The "solar Z" , "Chabrier IMF" and "KA97" results are obtained using stellar age and mass estimates with a 
different metallicity, a different IMF, and a different stellar population synthesis model from the baseline analysis, respectively. The 
"Ay < 0.5" results are obtained with a more strict cut about the dust extinction of host galaxies. The "SN ±0.3 mag" results are 
obtained with the SN la light curve ±0.3 mag fainter/brighter than in the baseline analysis. All errors are in statistical la. 
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Fig. 9. The SN la DTD foitla.) estimated by different prescriptions from the baseUne analysis are shown by hnes. The data points 
are for the baseline analysis, which are the same as those in Fig. 7. The labels for the lines are the same as in Table 4, and see this 
table and main text for explanations. 



prescriptions of the photometric redshift calculation as fol- 
lows. We first try the two cases of changing metallicity 
into the solar abundance {Z = 0.02) and changing IMF 
into the Chabrier (2003) IMF, from the baseline analysis. 
The choice of the supersolar metallicity {Z = 0.05) in the 
baseline analysis is motivated by the metallicity estimates 
of local early-type galaxies (Barber et al. 2007; Jimenez ct 
al. 2007), as mentioned in §2.1. Recent observations of the 
mass-metallicity relation at various redshifts give a further 
support to this choice. According to the mass-metallicity 
relation of Savaglio et al. (2005) measured for galaxies 
at a similar redshift range (0.4 < z < 1.0) to our sample, 
metallicity of the old galaxies selected in this work is ex- 
pected to be larger than the solar value, even when our 
definition of the stellar mass (integration of SFR) is taken 
into account. Therefore our examination in the metallic- 
ity range oi Z = (1-2.5) Zq is reasonable to check the 
systematic uncertainties in DTD. 

Furthermore, we used a completely different popula- 
tion synthesis model by Kodama & Arimoto (1997, here- 
after KA97), based on an independent stellar spectrum li- 
brary. We used the My = —21.98 model in the metallicity- 
sequence KA97 models as a template. In this model, a gas 
infall and galactic wind are taken into account rather than 
a simple exponential SFR evolution, and the chemical evo- 
lution is solved in a self consistent manner. This is a model 
for a local elliptical galaxy, and the star formation time 
scale is about tsf 0.1 Gyr. We have only one template 
for the star formation history, and hence the age estimate 



is likely to be less accurate than the baseline analysis, but 
we can check the dependence on different stellar popula- 
tion synthesis models as well as on the number of used SF 
history templates. 

In Fig. 10, we show the distribution of (t*}/(i*)bi and 
M,/M,. bi of the old galaxies, i.e., the ratios of (t,) and 
Ms, estimated by different prescriptions to those in our 
baseline analysis. The logarithmic means and standard 
deviations of these quantities are tabulated in Table 5. It 
can be seen that the means of are not significantly 
changed, and the typical dispersion of (<*)/(t*)bi is about 
a factor of two. This means that the bin width of our 
DTD measurements is not unreasonably small compared 
with the systematic uncertainties. As for stellar masses, 
systematic offsets of mean values of Af*/A/,.bi from unity 
can be seen. However, for the solar Z and the Chabrier 
IMF cases, these offsets are considerably reduced when 
we consider the K-hand luminosity at the age of 11 Gyr 
(Ljifi) as the mass indicator. The mass-to-light ratios are 
M^/Lk,o = 2.0 and 1.4 [Mq/Lk,q] for the solar Z and 
the Chabrier IMF cases, respectively (1.8 for the baseline 
analysis), and hence the effective offsets in log^oiif.o now 
become 0.056 and —0.043, respectively, i.e., less than 
10% difference in DTD estimates. This demonstrates the 
merit of using Lk,q as the mass indicator to reduce the 
uncertainties about stellar mass estimates. On the other 
hand, M^/Lka = 2.1 [Mq/L^,©] for the KA97 model and 
it cannot compensate the large offset of A/»/A/*.bi- This 
is most likely a result of the limited SF history template 
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in the estimates using the KA97 model (only 1 template 
corresponding to tsf ~ 0.1 Gyr). It is expected that the 
age estimates with a single small value of tsf tend to be 
underestimate of the true age, and this trend is in fact 
seen in the left panel of Fig. 10. The underestimate in 
stellar age would result in an underestimate in the stellar 
mass. 

To examine the sensitivity of our DTD estimates to 
these uncertainties of stellar ages and masses, we repeated 
the DTD estimates but using M*, {M^,/ Lka), and (tia) 
calculated with the different prescriptions, and the results 
are shown in Table 4 and Fig. 9. The change in the DTD 
is within the statistical Icr errors of the baseline results 
for most of the data points. The KA97 model case shows 
a large deviation from the baseline analysis at tia_ ^ 0.5 
Gyr, compared with the other cases. This is probably 
due to the age underestimates by a single SF template as 
discussed above, and we consider that the baseline results 
are more reliable than those using the KA97 model. The 
data points at tia = 0.5-2 Gyr are quite robust against 
these tests. 

Furthermore, we calculate (t*)spcc and M^^gpoc for the 
old galaxies with spectroscopic redshifts, to examine the 
systematic uncertainties in stellar age and mass coming 
from the uncertainty about redshift. These are calculated 
from the same photometric redshift calculations, but fix- 
ing the redshift at the value known by spectroscopic obser- 
vations. The results are shown in Fig. 11, as the ratios to 
(i*)phot and M^^phot obtained by the normal photometric 
redshift calculations (i.e., redshift treated as a free param- 
eter). The logarithmic mean and the standard deviation 
of these ratios are shown in Table 5. Many galaxies have 
exactly the same value of (t*)spoc = (i*)phot on the age 
grids of the hyperz code, and 73% of the old galaxies have 
(^*}spcc within a factor of two from (t*)phot- It seems that 
these uncertainties are not greater than the others about 
the photometric age/mass estimates discussed above. 

A caveat in our analysis is that the variety of SF his- 
tory used in photo-z calculation is limited. It is computa- 
tionally impractical to increase the number of SF history 
templates by a large factor from that used in this work 
(eight). However, the dependence on different SF histo- 
ries has partially been tested by the KA97 model case as 
described above. We also note that, since we selected old 
galaxies, the estimates of mean stellar age are expected 
to be rather insensitive to assumed SF histories. What 
is important for DTD estimates is the mean stellar ages 
rather than the details of SF history, provided that the 
dispersion of stellar age around the mean is not greater 
than the bin width of DTD measurements. This has al- 
ready been checked by the calculations of (trtia ) in section 
4.2. 

A considerable fraction of the old galaxies have Ay > 0.5 
in spite of low star formation activity, as shown in Fig. 1. 
It might be an artifact due to the uncertainties in Ay es- 
timates, but in that case it would affect the age estimates. 
Therefore we removed the old galaxies with Ay > 0.5 and 
repeated the DTD calculations, which are shown in Table 
4 and Fig. 9. The number of the SN candidates is re- 



duced to 48, but the change from the baseline analysis is 
not significant. The effect of dust extinction on supernova 
luminosities will be discussed in the next subsection. 

We have removed the nuclear region of host galaxies 
within 1.9 pix radius in our SN candidate search, and it 
might induce some biases in the age estimates if the re- 
moved nuclear regions have considerably different stellar 
age from the outer regions. The mean fraction of i' flux 
within the nuclear region of the old galaxies is 35% of 
the total isophotal flux, and hence the outer regions are 
more dominant than the nuclear regions in the photomet- 
ric redshift calculation. Figure 3 indicates that most of 
the old galaxies have elliptical morphologies, and they 
are likely to be the ancestors of the present-day ellip- 
tical galaxies (Yamada et al. 2005). It is known that 
the radial color gradient of local elliptical galaxies is ex- 
plained by metallicity gradient rather than age gradient 
(Tamura et al. 2000). The typical metallicity gradient 
is d\ogZ/d\ogr 0.3 (Kobayashi & Arimoto 1999), and 
the effect of metallicity gradient has already been tested, 
at least partially, by the above examination of the age 
dependence on metallicity. 

4-4- Systematic Uncertainties in Supernova Luminosity 

The light curves of SNc are essential for calculations of 
the visibility time and hence for a rate study. We used the 
standard light curves of various supernova types as in Oda 
& Totani (2005). However, we selected SNe la in galaxies 
with old stellar population, and a systematic trend that 
SNe la in old stellar population are fainter than those 
in star forming galaxies by a peak magnitude difference 
of AAIb ^ 0.3 has been known (Gallagher et al. 2005; 
Sullivan et al. 2006). 

Extinction in host galaxies may also change the appar- 
ent brightness of SNe, and hence the effective luminos- 
ity. We assumed a low extinction of E{B — V) = 0.05 (or 
Ay = 0.155 for the standard Milky Way extinction curve) 
for SNe la, which is reasonable for old galaxies. However, 
this value is smaller than (Ay) — 0.31 of the old galaxies 
estimated by the photo-z code, though this value should 
not naively be taken as the extinction of SN flux, as men- 
tioned in section 3.1. It should also be noted that the dis- 
tribution of old galaxy Ay is strongly peaked at Ay ~ 
(Fig. 1), and the uncertainty in Ay estimate may have 
resulted in a larger {Ay) than the real value. 

To examine the sensitivity of the DTD estimates to the 
effective changes of SN luminosity by these effects, we 
repeated the DTD estimates with the SN la light curve 
shifted by ±0.3 mag, and the results are shown in Table 
4 and Fig. 9. The difference from the baseline analysis 
is small, and the uncertainty about the SN luminosity 
or host galaxy extinction is unlikely to change our main 
conclusions. 

4-5. Uncertainties in the Comparison with the SN la 
Rate of Local Elliptical Galaxies 

We assumed 11 Gyr for the mean age of the local el- 
liptical galaxies and hence for the data of Mannucci et 
al. (2005). We have chosen this value as a mean value 
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Fig. 10. The histograms of the mean stellar ages (tt) (left) and stellar mass M» (right) of the old galaxies estimated by alternative 
photometric redshift calculations to the baseline analysis. They are shown in the ratios to those in the baseline analysis. The solid 
and dashed lines are for the cases of changing to the solar metallicity {Z = 0.02) and to the Chabrier IMF, respectively. The dotted 
line is for the estimates based on a completely different model (KA97) of stellar population synthesis. 



Table 5. The logarithmic mean and standard deviation of the ratios of mean stellar age and stellar mass to those in the baseline 
analysis, for alternative photometric redshift calculations using the solar metallicity, the Chabrier IMF, the KA97 galaxy evolution 
model, and fixing redshifts at the spectroscopic values. 
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*These quantities should be replaced by (t»)spcc/(t»)phot ^i^d -'^^*,spcc/A/, p^ot for the "spec" column. 




Fig. 11. The "spectroscopic" mean stellar ages ((t,)) and stellar masses (M,) versus spectroscopic redshifts, for the old galaxies 
with available spectroscopic data. See the main text for the definition of these quantities. They are given as the ratios to those by 
the normal photometric redshift calculations. 



16 



Totani et al. 



of the luminosity-weighted stellar ages of the early type 
galaxies at 0.05 < z < 0.2 estimated by Jimenez et al. 
(2007), based on the MOPED fits to the SDSS spectra. A 
slightly smaller but similar age (~ 9 Gyr) is obtained for 
the SDSS luminous red galaxies at 0.15 < z < 0.4 (Barber 
et al. 2007), when the look back time to the redshift of 
these galaxies is corrected. The age dispersion around the 
mean for these early-type galaxies is about 2 Gyr, both 
for the samples of Jimenez et al. and Barber et al. A sys- 
tematic uncertainty of 2 Gyr in the 11 Gyr data point 
does not significantly affect the conclusion that the DTD 
is well described by a power-law. ^ 

In the SN la rate estimate for the local elliptical galax- 
ies (Mannucci et al. 2005), very faint SNe la with the 
stretch factor s < 0.8 are included, reflecting the observed 
trend that SNe la found in old galaxies are subluminous 
on average (Gallagher et al. 2005; Sullivan et al. 2006). 
On the other hand, such subluminous SNe la are hardly 
detected in flux-limited high-^ supernova searches [see, 
e.g., Fig. 11 of Sullivan et al. (2006)]. If the fraction of 
subluminous SNe la does not evolve with the delay time, 
we should correct our DTD data points upward to make 
an unbiased comparison with the data of Mannucci et al. 
(2005). However, the fraction of s < 0.8 SNe la in local 
elliptical galaxies is ~ 30% [see again Fig. 11 of Sullivan 
et al. (2006)], and it does not significantly affect the main 
conclusions of this work. Furthermore, since our SN can- 
didate sample is at high redshift, it mainly probes shorter 
delay times of tia ^ 2 Gyr, and there may not be sublu- 
minous SNe la with such short delay times. In this case, 
the correction discussed here is not necessary. 

5. Implications for the SN la Progenitor 

5.1. the Double- Degenerate Scenario 

The DTDs predicted by the four different theoretical 
models (Ruiz-Lapuente & Canal 1998; Yungelson & Livio 
2000; Greggio 2005; Belczynski et al. 2005) based on the 
DD scenario are shown in Fig. 12 in comparison with 
the observed DTD data points. The model curves are 

^ Very recently, Gallagher et al. (2008) reported the age and 
metallicity estimates for the elliptical host galaxies of the lo- 
cal SNe la, and the ages estimated by them have a consider- 
ably wider distribution than those of Jimenez et al. (2007) and 
Barber et al (2007), including a significant fraction of young 
galaxies with ages less than 5 Gyr. There is no difference in the 
properties of the SN la host galaxies and field elliptical galaxies 
(Gallagher et al. 2008), indicating that the age difference comes 
from some systematic uncertainties in the age estimates. The 
origin of the difference is not clear to us, but it should be noted 
that the age estimate by Gallagher et al. (2008) is based on H/3, 
FcA5270 and FcA4383 indices with the spectral library including 
only single-burst stellar population (SSP). In such analyses, the 
age estimates are easily affected by a small amount of recent 
star formation, and they stated in their paper that their age 
estimates should be regarded as lower limits to the true ages. 
Their figure 7 shows a strong correlation between the ages and 
metallicities especially for galaxies younger than 5 Gyr, indicat- 
ing an effect of the age-metallicity degeneracy in the fittings. 
The estimates by Jimenez et al. (2007) are based on full spec- 
tral information including continuum as well as absorption line 
indices, allowing any SF history in If time bins. 



normalized at the data point of tja = 11 Gyr, because the 
DTD behavior at large delay time is especially important 
for this scenario from the theoretical point of view (see 
below). However, the normalization should be regarded 
as a free parameter and we should compare only DTD 
shapes between the observations and the models. It is 
impressive that the predictions by different authors are 
very close, and in excellent agreement with the observed 
DTD at iia ^ 0.2 Gyr. It is not surprising that different 
authors made similar predictions, because it is a robust 
and generic prediction of the DD scenario as argued below. 

A common feature for DTD models based on this sce- 
nario is that the delay time is mainly determined by the 
time tow from the formation of a DD binary (i.e., a bi- 
nary of two white dwarfs) to a merger after the angular 
momentum loss by gravitational wave radiation, especially 
in the long delay time range of tia ^ 1 Gyr. The general 
relativity tells us tja ~ ^gw oc a^, where a is the initial sep- 
aration of the DD binary. If the separation distribution is 
given by /scp(a) oc a^, the DTD should be 

foih.) cx /3ep(a) cx ti}'-^^/'' , (7) 

atia 

and hence a = —(3 — /3)/4. It is known that the distri- 
bution of initial binary separation aq at the time of bi- 
nary star formation is approximately flat in logag, i.e., 
/scp,o(ao) Of 0,0^ (Abt 1983). Although there is a signif- 
icant change and contraction of binary separations dur- 
ing the evolution from a binary star formation to a DD 
binary (e.g., Greggio 2005), it would be reasonable to as- 
sume P —1 also for the initial separation of DD binaries. 
Then we expect a ~ — 1, as found in the DTD models 
based on the DD scenario. Though (3 is rather uncer- 
tain and model-dependent, the dependence of a on /3 is 
small. Furthermore, the range of the delay time from 0.1 
to 10 Gyr corresponds to a range of a only by a factor of 
3.2. Therefore, if a is smoothly distributed in this nar- 
row range, we expect a power-law DTD with a — 1 in 
wide and general conditions. This is a general result ap- 
plicable to any merging phenomena triggered by gravita- 
tional wave radiation, such as binary neutron star mergers 
(Totani 1997), and this is why the DTD models by differ- 
ent authors are quite similar to each other, especially at 
large delay times. The agreement of the measured DTD 
with this generic prediction gives a strong support for the 
DD scenario. 

In Fig. 12, we selected the standard DTD model 
when a few models based on the DD scenario are pre- 
sented by a single author group [the standard DDS model 
of Belczynski et al. (2005) and the close-DD model of 
Greggio (2005)]. The variation models predict slightly 
different value of a, but smooth power-law like DTDs are 
always common predictions, as expected from the above 
consideration. 

On the other hand, the difference between the DTD pre- 
dictions by different authors becomes bigger at short delay 
time (tia ^0.1 Gyr), where the delay time is dominated 
by stellar evolution time scale rather than <gw- There are 
large uncertainties in the treatments of stellar evolution in 
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Fig. 12. The observed SN la DTD foitis.) compared with the theoretical predictions based on the DD scenario by Ruiz-Lapuente & 
Canal (1998, dotted), Yungclson & Livio (2000, dot-dashed), Grcggio (2005, solid), and Belczynski ct al. (2005, dashed). The data 
points are the same as those in Fig. 7. The model curves are normalized by the DTD data at tja = 11 Gyr, but the normalization 
is arbitrary and only the DTD shapes should be compared between the data and models. 



binary population synthesis, and the discrepancy between 
some of the DTD models and the observed data at this 
short range is not serious for the DD scenario. 

It is generally assumed that DD binaries are formed 
through the common envelope evolution phase. If, in addi- 
tion, DD binaries are formed also from intermediate mass 
binaries through the classical Roche lobe outflow phase, 
DD binaries of this channel with large separations may 
affect the late time behavior of the DTD in the DD sce- 
nario (De Donder & Vanbeveren 2003). In other words, 
the observed DTD gives some constraints on the efficiency 
of DD binary formation via the Roche lobe outflow phase. 

5.2. the Single-Degenerate Scenario 

The measured DTD is compared with the theoretical 
models based on the SD scenario in Fig. 13. Again, the 
normalization of the models is arbitrary, and we should 
compare only the shape of the DTD function. Here, the 
theoretical curves are normalized by minimization to 
the data (our own measurements at 0.1-8 Gyr plus the 11 
Gyr data point). In contrast to the DD scenario, the pre- 
dictions by this scenario are quite different depending on 
different authors. It is common to this scenario that the 
delay time is essentially determined by the main-sequence 
life time of the secondary star in a binary, because other 
time scales (e.g., accretion phase onto the white dwarf) 
are much shorter. However, the DTD models shown in 
this figure arc calculated by considerably different meth- 



ods with a wide range of complexity. 

The predictions by Matteucci et al. (2006) and the 
Z = 0.05 model of Kobayashi & Nomoto (2008, hereafter 
KN08) are close to a simple power-law, and they arc rel- 
atively in good agreement with the observed DTD, com- 
pared with other models. However, it should be noted 
that the two models are based on the simplest calcula- 
tions among the models shown here, though such a simple 
analytical approach is useful to see the general behaviors 
and in applications for e.g., studies of chemical evolution 
in galaxies (Kobayashi ct al. 1998; Matteucci & Rccchi 
2001; Matteucci et al. 2006; KN08). In these calcula- 
tions, the condition for a binary system to evolve to a SN 
la event (i.e., stable accretion up to the Chandrasekhar 
mass) is determined by the initial masses of the primary 
and secondary stars, regardless of the binary separation. 
If the masses are in acceptable ranges, a binary evolves 
into a SN la with a constant efficiency, and no SN la is 
produced in other mass ranges. In such a calculation, the 
DTD shape is determined by IMF and the distribution of 
the primary /secondary mass ratio. If smooth functions 
are assumed for these, a smooth DTD without a charac- 
teristic scale of tia is obtained. However, a constant SN la 
efficiency within an extended mass range is likely to be an 
oversimplification. The SN la condition should be deter- 
mined by the combination of the stellar masses and binary 
separation, and it is natural to expect that there arc some 
characteristic secondary mass scales preferred for success- 



18 



Totani et al. 




Fig. 13. The same as Fig. 12, but for a comparison with the DTD predictions based on the SD scenario. The red curves are 
predictions of Greggio (2005, solid curve), Matteucci et al. (2006, dotted), and Kobayashi & Nomoto (2008, dashed), based on 
analytic calculations in which the condition for successful SNc la is determined by stellar masses in a binary. The blue curves are 
predictions of Ruiz-Lapucnte & Canal (1998, dotted), Yungclson & Livio (2000, dot-dashed), Belczynski et al. (2005, dashed), and 
Meng et al. (2008, solid), based on binary stellar population synthesis calculations. The theoretical curves arc normalized by 
minimization to the data, and only the DTD functional shapes should be compared in this plot. 



ful SN la events, which should appear as particular time 
scales in DTD, rather than a power-law like DTD. 

In the KNOB model based on the progenitor model of 
Hachisu et al. (1996, 2008a), there are two distinct pop- 
ulations of the main-sequence donor channel (WD-I-MS) 
and the red-giant donor channel (WD-I-RG) for short and 
long delay times, respectively. The relative abundances of 
the two are determined empirically by fitting to the chem- 
ical evolution data, rather than by theoretical modeling. 
Therefore, the good agreement of the KN08 model with 
our DTD estimate may partially be a result of this feed- 
back from other observations, indicating that our DTD es- 
timate is also consistent with the chemical evolution data. 
Though the WD-I-MS channel is widely considered as the 
promising SD progenitor of SNe la, some theoretical cal- 
culations indicate that the SN la rate by the WD-I-RG 
channel is much lower than that of the WD-I-MS chan- 
nel (Yungelson & Livio 1998; Han & Podsiadlowski 2004; 
Meng et al. 2008, but see also Hachisu et al. 2008b). 

The SN la condition is determined only by stellar 
masses also in the standard SD-Chandra model of Greggio 
(2005), but the mass budget of accretion up to the 
Chandrasekhar mass is treated in more detail. An in- 
teresting feature is the sharp drop of DTD at a large 
delay time of 8 Gyr. SNe la with such a large delay 
time are produced from binary systems with a small sec- 
ondary mass. Such a low mass star has a small enve- 



lope mass that can be used for accretion onto its compan- 
ion. Therefore, the primary stellar mass must be large 
enough to ensure that the initial white dwarf mass is suf- 
ficiently massive and the white dwarf successfully grows 
to the Chandrasekhar mass. Consequently, the number of 
binaries that can evolve to SNe la rapidly decreases with 
decreasing secondary mass. It seems difficult to reproduce 
a single power-law like DTD up to '--^lO Gyrs if this effect 
is incorporated. 

Finally, predictions of the SG-Ch model of Yungelson 
& Livio (2000), the standard SDS model of Belczynski 
et al. (2005), the CLS model of Ruiz-Lapuente & Canal 
(1998), and the model of Meng et al. (2008) (acE = 3 
and Z ~ 0.03) are shown, which are based on detailed cal- 
culations of the binary population synthesis, where the 
calculations start from the initial conditions including the 
primary and secondary stellar masses and the binary sep- 
aration. There are a number of parameters and uncertain 
physical processes in such calculations, and different con- 
ditions for SNe la are assumed by different authors, lead- 
ing to vastly different DTD predictions. However, a clear 
trend is that some characteristic scales of the delay time 
appear in the DTD. making the DTD shape more complex 
than the simple analytic models. This is reasonable in re- 
alistic conditions, as argued above, and it is theoretically 
unlikely that a simple power-law like DTD is obtained in 
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the framework of the SD scenario^. 

6. Conclusions 

6.1. Summary of This Work 

We measured the delay time distribution (DTD) of type 
la supernovae by using the statistics of the faint variable 
objects detected in the systematic variable object survey 
performed as a part of the Subaru/XMM-Newton Deep 
Survey project. 

Based on photometric redshift calculations using 9 band 
{BVRci'z'JK, 3.6 /im, and 4.5 ^m) photometries, we se- 
lected 16,492 old galaxies from the SXDS data in the red- 
shift range of 0.4 < z < 1.2, by requiring that their SED can 
be fit by a single starburst and their ages are significantly 
greater than the star formation time scales of the bursts. 
This selection is essential to our work, because stellar age 
is expected to be a good estimator of la delay time in 
such galaxies. Furthermore, we have two more merits of 
this selection: (1) delayed SNe la are dominant compared 
with CC SNe and prompt SNe la, and (2) extinction effect 
should be small. 

We then selected variable objects associated with these 
galaxies but having significant offsets from the nuclei of 
galaxies, to remove the contamination of AGNs. We found 
65 variable objects, whose locations with respect to their 
host galaxies closely trace the profiles of galactic light, and 
hence the majority of them must be supernovae. Though 
we do not have spectroscopic confirmation of the SN types, 
we quantitatively demonstrated that the majority (^80 
%) of the SN candidates should be SNe la, based on the 
variability luminosity, redshift distribution, and proper- 
ties of the host galaxies. 

Then the DTD in the delay time range of tia = 0.1- 
8 Gyr is derived by calculating the delay time for each 
SN la candidate from age and star formation history of 
host galaxies. Combined with the observed SN la rate 
in elliptical galaxies in the local universe, we derive the 
SN la DTD in a range of tia = 0.1-11 Gyr, and found 
that it can well be described by a simple power-law, 
fDiti^)=0-55t°o.ll{tiJl Gyr)« century-i(lOi°LK,o,0)-' 
with a = -1.08^0,11. Here, the DTD function /n(tia) is 

^ After wc put the first version of this paper on the preprint server 
(arXiv:0804.0909vl), Hachisu, Kato, & Nomoto {2008b) per- 
formed a DTD calculation based on binary population synthesis 
incorporating their progenitor model. An important difference 
of this calculation from KN08 is that the ratio of the two evolu- 
tionary channels (WD-I-MS and WD+RG) is calculated theoret- 
ically rather than determined by fitting to observed data. They 
found that their DTD prediction is in good agreement with the 
data derived here, indicating that there is a viable model param- 
eter space in the SD scenario. On the other hand, each of the two 
channels has a narrow DTD shape similar to the other SD mod- 
els based on binary population synthesis, confirming the trend 
of the SD predictions. It seems somewhat a fine-tuning that two 
independent components form a featureless power-law DTD as 
a sum, although it could well happen in the nature. The ratio 
between the two components depends on still uncertain physi- 
cal processes, such as mass-stripping of secondary stars by wind 
from accreting white dwarfs. Further theoretical investigation 
is important. 



per unit delay time and per unit mass of a single-burst 
stellar population, and Lk.o (A'-band luminosity at an 
age of 11 Gyr) is used as an observational estimator of 
the stellar mass. We performed various tests about the 
systematic uncertainties in this DTD measurement, but 
the changes of DTD estimates are not large enough to 
change our main conclusions significantly. We tried a va- 
riety of DTD calculations with different prescriptions as 
shown in Fig. 9, but the decreasing trend of /^(tia) con- 
sistent with foitia.) at tia ^ 0.3-3 Gyr is found in 
all cases, and the data points at 0.5-2 Gyr are especially 
robust against examined systematic uncertainties. 

The derived DTD at ti^, ^ 0.2 Gyr is in excellent agree- 
ment with the theoretical predictions based on the DD 
scenario. The theoretical predictions by different authors 
are very similar to each other, and a featureless power-law 
shape is inevitable consequence of the general relativity in 
this scenario. Therefore we consider that the agreement 
between the observed and predicted DTDs gives a strong 
support to this scenario. It indicates that the major con- 
tributor to SNe la is the DD channel for delay times larger 
than tia ^ 0.1 Gyr, although some contribution from other 
channels cannot be excluded. 

On the other hand, the predictions by the compet- 
ing SD scenario are vastly different for different authors. 
Although the predicted DTDs based on simple analytical 
approaches have smooth shapes which are broadly consis- 
tent with our measurement, predictions based on more 
detailed binary population synthesis calculations show 
strong peaks in the DTD shape, which do not fit the ob- 
served DTD. This trend can naturally be understood; if 
there is any preferred scale of the secondary stellar mass 
for a binary to successfully evolve into a SN la, it should 
be reflected as a characteristic delay time scale in this 
scenario. Therefore we consider that our result does not 
favor the SD scenario in general as the major channel to 
the delayed SN la population (ija ^0.1 Gyr). However, 
there are many degrees of freedom and uncertainties in 
the theoretical modeling based on the SD scenario, and it 
would be premature to reject the SD scenario simply from 
our result. Our result should set strong constraints on the 
model parameter space if the SD channel is responsible for 
the majority of the delayed SNe la, and it would be use- 
ful for the future theoretical studies of the SD progenitor 
models. It is highly desirable to examine physical effects 
and evolutionary paths that are not taken into account in 
existing models of binary population synthesis. 

6.2. Discussion and Other Implications 

In the literature, it has occasionally been argued that 
the SD scenario is more favored than the DD scenario, 
but arguments against the DD scenario are not particu- 
larly strong. A merger of two white dwarfs may result in 
an accretion induced collapse rather than a SN la (Saio 
& Nomoto 1985), but theoretical uncertainties are still 
quite large (Piersanti et al. 2003; Yoon et al. 2007). 
Many binary systems are proposed as candidates of the 
SN la progenitor in the SD framework, compared with 
the observed number of DD binaries having total masses 
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larger than the Chandrasekhar mass (Parthasarathy et al. 
2007). However, DD binaries are difficult to detect and 
the statistics of DD binary searches is not sufficient yet 
to confirm or reject the DD scenario (Tovmassian et al. 
2004; Geier et al. 2007; Napiwotzki 2007). There are 
some implications for the progenitor from studies of SN 
la remnants (Ruiz-Lapuente et al. 2004; Badenes et al. 
2007; Ihara et al. 2007), or from spectroscopic studies 
(Leonard 2007; Patat et al. 2007a, b; Simon et al. 2007) 
and archival progenitor searches (Voss & Nelemans 2008; 
Nelemans et al. 2008; Roelofs et al. 2008) for nearby 
SNe la. However, conclusive results from these methods 
have not yet been obtained, mainly because of the limited 
statistics and theoretical uncertainties. 

Another possibility is the SD sub-Chandra scenario 
(also known as helium ignitors or edge- lit detonations), in 
which white dwarfs explode as SNe la before they reach 
the Chandrasekhar mass. Since the required amount of 
accreting mass is smaller than in the SD Chandra sce- 
nario, binaries with lower secondary masses can more eas- 
ily evolve to SNe la. This might be useful to overcome 
the difficulty of the SD scenario to explain SN la rate at 
tia ^ 10 Gyr (Greggio 2005). However, this scenario is cur- 
rently not popular because predicted spectra are in serious 
disagreement with observations (Hillebrandt & Nicmcyer 
2000; Livio 2001). 

Recent studies (Scannapieco & Bildsten 2005; Sullivan 
et al. 2006; Aubourg et al. 2007) have tried to model the 
DTD by two components: the prompt component propor- 
tional to SFR and the delayed component proportional to 
stellar mass (i.e., constant /d). Although our result does 
not put a strong constraint on the amount of the prompt 
component, a constant fr) seems to be an oversimplifica- 
tion for the delayed component, since foitia) at an inter- 
mediate delay time of tja ~ 1 Gyr is about 10 times larger 
than /d(11 Gyr). It should be noted that the prompt la 
population inferred from our data is a considerable frac- 
tion (~ 20%) of all SNe la integrated over ti^ = 0-11 Gyr, 
if we make a modest assumption that foitia.) at tja < 0.1 
Gyr is constant at the value of /£)(0.1 Gyr). Most of the 
observed SN la rate dependence on galaxy properties can 
be reproduced by a DTD similar to those predicted by the 
DD scenario (Greggio 2005; Mannucci et al. 2006). The 
observed enhancement of SN la rate in radio galaxies may 
require an even higher prompt la fraction than that ex- 
pected from our data or the DD scenario (Delia Valle et al. 
2005; Mannucci et al. 2006), if the enhancement is due to 
a recent starburst connected to AGN activity. However, 
this interpretation seems to be inconsistent with the fact 
that no enhancement of CC SN rate is observed in radio 
galaxies (Greggio, Renzini, & Daddi 2008). The statistics 
of this enhancement is still small {2a level), and it must 
be confirmed by future observations. 

Recently, Pritchet et al. (2008) reported that DTD pro- 
portional to ij"" '^^"'^ is implied based on the SNLS data 
(Sullivan et al. 2006), which seems to be inconsistent with 
our results. It should be noted that their constraint on 
DTD is not based on delay time estimate for each SN la, 
but it is indirectly derived from the correlation between 



SN la rate per unit galaxy mass (specific SN la rate) and 
SFR per unit galaxy mass (specific SFR). It is difficult to 
estimate how large is the systematic uncertainty in such 
an analysis. Especially, low specific SFR galaxies should 
be treated with caution. There should be no tight rela- 
tion between specific SN la rate versus specific SFR in 
such galaxies; specific SFR could change significantly by 
changing SFR with a fixed stellar mass, but specific SN 
la rate is hardly affected if it is dominated by old stellar 
population. This should limit the power of this approach 
to constrain DTD. It should also be noted that the stel- 
lar mass estimates in Sullivan et al. (2006) are based 
on five optical bands {u*g'r'i' z'), although near-infrared 
bands are essential to reliably estimate stellar masses of 
galaxies. Another important point is that we obtained a 
strong constraint on the slope index a by supportively us- 
ing the SN la rate in nearby elliptical galaxies. If we fit 
only our own DTD data, we obtain a — — 0.921q27. It is 
important to combine the nearby rate data (large tia) and 
high-z data (small <ia) to get a strong constraint on the 
DTD shape in a wide range of the la delay time. 
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